Vorticity Generation in Slow Cooling Flows 
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We show that any generic non-adiabatic slow flow of ideal compressible fluid develops a significant 
vorticity. As an example, an initially irrotational conductive cooling flow (CF) is considered. A 
perturbation theory for the vorticity generation is developed that employs, as a zero order solution, 
a novel two-dimensional similarity solution. Full gasdynamic simulations of this CF demonstrate 
the vorticity generation and support the theory. The relevance of this problem to the experiments 
with the "hot channels" is discussed. 
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The most general mechanism of the vorticity produc- 
tion in non-adiabatic flows of ideal compressible fluids 
relies on the misalignment of pressure and density gradi- 
ents H. Recent experiments with supersonic flows [pP] 
have clearly demonstrated the efficiency of this mecha- 
nism (which is called baroclinic). The baroclinic mecha- 
nism can also operate in slow gas flows, and its speciflc 
manifestations in meteorology m and combustion B are 
known. The main objective of this Letter is to show that 
the vorticity production represents a generic and signifi- 
cant property of any slow non-adiabatic gas flow. 

For concreteness, we will consider the conductive cool- 
ing flows (CFs) and refer to the "hot channels" produced 
in the air by lasers or electric discharges |g,||. After 
pressure equilibration these channels develop a signifi- 
cant vorticity and small scale turbulence and cool much 
faster than they would because of molecular thermal con- 
duction. Picone and Boris [Q interpreted these results in 
terms of the baroclinic vorticity production during the 
rapid channel expansion (that is, on the acoustic time 
scale) M . Schlicrcn photographs of the hot channels pP] 
clearly show that the most significant vorticity dynamics 
occurs on a much longer time scale. According to Picone 
and Boris, "after pressure equilibration... vorticity is no 
longer generated, however, significant residual vorticity 
exists" 0. We wish to present an alternative scenario 
which assumes that a significant vorticity is created on 
the long, heat-conduction time scale. The underlying 
physics is the following. After a few acoustic times, fol- 
lowing the rapid energy release, the gas pressure becomes 
very close to the (constant) ambient pressure, while the 
vorticity generated earlier is presumably damped out. As 
the temperature inside the channel is still very high, a 
low-Mach-number conductive CF develops that cools the 
channel by filling it with the cold gas from the periphery. 
Slow conductive CFs were studied previously in the con- 
text of a "point-like" energy release, like a high-altitude 
explosion in the atmosphere |9| or "fireball" produced by 
a laser spark in front of condensed matter jl^l ■ We will 
show that, unless the energy release geometry is fully 
symmetrical, small pressure gradients, intrinsic in the 
CF, result in a significant vorticity production. 



The simplest conductive CF of a perfect gas is de- 
scribed by the standard gasdynamic equations pT[|: 
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where d/dt = d/dt -I- v • V is the total derivative. The 
distance is measured in the units of a characteristic spa- 
tial scale of the problem tq (see later), while the time 
is measured in the units of the heat conduction time 
To = 7(7 — l)~^-RgPo^o'*~^^o~^- Furthermore, the gas 
density p and temperature T are scaled by their (pre- 
sumably constant) values "at infinity" po ^^nd Tq, the 
velocity v is scaled by ro/ro, and the pressure p = pT is 
scaled by {Rg/p)poTo. The non-dimensional parameter 
e — ro/cgTo (where c^ — RgTo/p) represents the char- 
acteristic Mach number of the flow. Finally, i?g,7 and 
p are the gas constant, adiabatic index and molar mass, 
respectively, while the heat conductivity is assumed to 
be a power-like function of the temperature: kT'^ in the 
scaled units, k = const. (For the molecular air v = 1/2.) 

We start with a perturbation theory that describes the 
initial stage of the vorticity production. Then we report 
on numerical simulations with the full equations (|l|)-(|^) 
that show the vorticity generation in the same CF and 
support the theory. 

In the low Mach number regime, e^ ^ 1, the tem- 
perature and density contrasts can still be large, but 
pressure non-uniformities are already small: p(r, t) — 
1 + e^ Sp{r, t) . Then, neglecting the small 5p terms in Eq. 
(^) and equation of state, we obtain V • (v — T'^VT) = 
and pT — 1, respectively. It follows that v — Vp -I- Vg, 
where Vp = —p^'^^^'SI p is the irrotational component of 
the fluid velocity, and v., is the solcnoidal component: 
V • Vg =0. Substitution of v into Eq. (1^) yields a non- 
linear transport equation: 



^ + (vs-VV = V.(p-''-Vp). (4) 

An important additional equation follows from Eq. (0): 

where a; = Vxv = Vxvsis the vorticity. Eq. (||) is 
equivalent to the well-known vorticity equation |ll| , as its 
right-hand-side can be rewritten as Vp x VSp/p^. Note 
that Eq. (||) does not include e, therefore, the vorticity 
production rate is, in general, of order unity. 

In this Letter we address the vorticity production in 
an initially curl-free flow. Accordingly, we assume that 
Va ^ Vp and, in the zero order, neglect the second term in 
the left side of Eq. (0) . The resulting nonlinear diffusion 
equation. 
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describes such a curl-free CF completely 
consider a first-order version of Eq. (0) 
terms of the vector field a(r, t) = dvs/dt: 
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V X aH X a = 
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rewritten in 
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with p and Vp given by the curl- free solution |T2|. Again, 
Eq. (R) shows that for a generic CF the vorticity produc- 
tion rate is of order unity. Therefore, the solenoidal part 
of the velocity field finally becomes comparable to its ir- 
rotational part (at which stage this perturbation scheme 
breaks down). 

Let us concentrate on a two-dimensional (2d) flow in 
the xy-plane with no z-dependence, where one can pro- 
duce the flrst two "classes of asymmetries" of the hot 
channels 0: (i) off-center laser beam propagation and 
(ii) non-circular cross section of the beam. Introduce a 
modified stream function ip^x, y, t), so that Ux — —dip/dy 
and Gy = dip/dx. Eq. (|^) becomes a scalar equation for 



V • (pV^) = 



^ + (vp.V)vp 



X Vp 
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where e^ is the unit vector in the z-direction. 

One should, however, deal first with Eq. (||) and find 
the zero-order solutions p(x,y,t) and Vp(a;,y,t) enter- 
ing Eq. (||). Remarkably, Eq. (||) has a family of 2d- 
similarity solutions of the second kind: 



p{x,y,t) =t'^+!^ R{i,ii), 



(9) 



where ^ — x/v , rj ~ y/t^, and f3 is an arbitrary real pa- 
rameter [nsl. Selection of parameter /3 requires the use 
of initial or boundary conditions. We shall adopt the fol- 
lowing initial density profile: p{x, y, t = 0) = A^^r''f{(f>), 



where r and (/> are the polar coordinates in the plane xy, 
and A and k are constants. In the case of a cylindrically- 
symmetric [/(</>) — 1] explosion along the z-axis, this 
profile with k — 2/(7 — 1) represents the r —>■ den- 
sity asymptotics that sets in at the end of the expansion 
The function /(</>) describes asymmetry, 
we extend this initial condition 



stage 
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In analogy to Ref. 
to the whole CF region. This idealization is justified (see 
Refs. 0,D and gasdynamic simulations below) as long as 
the density (temperature) contrast in the system remains 
much larger than unity. 

The initial condition yields /3 = {vk + k + 2)~^. Using 
Eq. (g), one arrives at a nonlinear elliptic equation for 
the shape function R(^,r]): 
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(we got rid of the constant A by choosing tq — A'^/^). We 
assume for simplicity that the initial density profile [and, 
hence, i?(^, 77)] is symmetric with respect to each of the 
Cartesian axes and solve Eq. (|l^) in the first quadrant 
with the no-flux boundary conditions at the ^- and rj- 
axes. In addition, we must require that i?(^ — > 00,77 ^^ 
00) = f^ficf)), where f and (j) are the polar coordinates 
in the plane £,, 77. Fig. 1 shows R{£,, rj) found numerically 
in a finite square for f{<j)) = 1 + a cos 20. We took the 
usual values 7 = 1.4 and v = 0.5 for the molecular air, 
and chose a — 0.6. (In this case /3 ~ 2/19, while the gas 
density at the channel axis grows in time like i^°/^^.) 
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FIG. 1. Contours of \gR{S^,ri). On the inner line 

IgR = 0.25 and increases outwards in steps of 0.25. Also, 
R(0,0) ^ 1.99. 

Now we return to Eq. (H). A similarity solution for 
p implies a similarity solution for ijj, that is il){x^y,t) = 
i""*(C, V), where a = 2{vk + k + l)/{vk + k + 2). For 



the shape function ^(^,?7) one obtains a hnear elhptic 
equation V ■ (i?V*) = (W x Vi?) • e^, where 



W = -{:^k + k + 2y 



dV dV 
,„* + * + l)V + «- + „- 

+ (V-V)V, 



V = —R^'^^^VR, and the V-opcrator now involves dif- 
ferentiation with respect to ^ and 77. This equation should 
be solved in the first quadrant with the Dirichlet bound- 
ary condition. We solved it numerically using the shape 
function R found earlier (in this case a — 34/19 w 1.79). 
The result is shown in Fig. 2. Now we can evalu- 
ate the vorticity w = ^(x, y, t) Gz, using the relation 
{diji!/dt)x,y = V^V- Following Picone and Boris 0, we 
calculate the vorticity flux f2 through the first quadrant 
as a function of time. The growth rate of this quantity, 
dfl/dt, is equal to 



/•OO /"OO /> 

/ / dxdy V^V = f (Vi/- • n) dl, 
Jo Jo Jc 



(11) 



where C is the contour going from infinity to zero along 
the y-axis and continuing to infinity along the x-axis and 
n is the external normal. Employing the similarity solu- 
tion for t/), we obtain 



dfl 



-jT = -it + toy 



(9* 



(12) 



where we have used the invariance of the similarity solu- 
tion with respect to a time shift and introduced to, the 
only fitting parameter of the theory. Integrating Eq. jl^ ) 
with a zero initial condition, we arrive at 



\n{t)\ = \B\{a-l)-^[tl-''-{t + toy 



(13) 



where B is the constant given by the expression in the 
square brackets in Eq. (12). Eq. ( p^ predicts a lin- 
ear growth of Q with time followed by saturation at a 
constant value \ B \ {a — 1) 
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FIG. 2. Contours of *(C,??)- On the outer line * = 10"" 
and increases inwards in steps of 10~". 



Now we report on the 2d-numerical simulations with 
Eqs. (||)-(||)- We employed an extended version of 
the code VULCAN |1^ that uses flexible moving grids 
and can operate in any combination of Eulerian and 
Lagrangian modes. In the rezoning stage we used the 
scheme of Van Leer |1^ that preserves second order ac- 
curacy. The code could work in an implicit mode, thus 
eliminating the Courant-Friedrichs-Lewy restriction on 
the time step. The initial conditions were 



p{r, 0, t = 0) = Pi, 



(5 + r'=(l + acos20) 
1 + r'=(l + acos2(/)) 



(14) 



for the density, v(r, (j),t — 0) ^ ^P^n'^^ Pin for the ve- 
locity, and unity for the pressure. For (5 <IC 1, the ini- 
tial density profile has an extended part described by 
r^(\ -|-acos2(/)) (which yields the 2d-similarity solution). 
On the other hand, pi„ is non-zero at r = and ap- 
proaches unity at r — > 00 as it should. In most of simula- 
tions we took 8 — 10^'^, a = 0.6 and e in the range of 10^^ 
to 10"^. Simulations show that the density history at the 
channel axis is described very well by the similarity scal- 
ing 1.99 (< + 4.2 X 10-5)i°/i9 until the late stage, when the 
density contrast is reduced. However, the velocity field 
(that was curl-free in the beginning. Fig. 3a) develops 
a noticeable vorticity which spatial structure is similar 
to that shown in Fig. 2. Finally, a distinctive vortex, 
advected towards the origin by the overdense gas inflow, 
appears (Fig. 3b). Since the problem is symmetric with 
respect to each of the axes, the corresponding "full" flow 
develops four symmetric vortices. Fig. 4 shows the time 
history of the vorticity flux through the first quadrant, 
ri, as found from the simulations. It is seen that the vor- 
ticity reaches a significant value. One can also see that 
the perturbation theory [Eq. (p^)] underestimates the 
saturated vorticity flux. This is understandable, as the 
perturbation scheme fails at large times. Interestingly, 
the agreement improves for a smaller value of io- 

In summary, we claim that any generic non-adiabatic 
gas flow develops a significant vorticity. For the low- 
Mach-number conductive CF that we have considered in 
detail, the further vorticity dynamics (instability?) is ap- 
parently sensitive to geometry (like in the Picone-Boris 
scenario) . We did not observe turbulence or other signifi- 
cant modification of the bulk transport properties in this 
(still highly symmetric) 2d-fiow. Correspondingly, the 
"hot channel" riddle requires further investigation. One 
can expect turbulence to show up in a less symmetric 3d 
situation, when perturbations along the channel axis are 
introduced. 

B.M. acknowledges a valuable discussion with P.V. 
Sasorov. 
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Density and velocity fields at f = (a) and 
^ (b). The velocity field (arrows) is scaled by 



FIG. 3. 
i = 3 X 10 

3000 (a) and 100 (b). On the inner density isolines p 
and increases outwards in steps of 0.1. 
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FIG. 4. Vorticity flux through the first quadrant vs time, 
as predicted by the full simulations, and by Eq. (|l^) with 
different values of to ■ 
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